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Abstract 

A general model for stationary, time-wise turbulent velocity is presented and discussed. 
This approach, inspired by modeling ideas of [2], is coherent with the K41 hypothesis 
of local isotropy, and it allows us to separate second-order statistics from higher order 
ones. The model can be motivated by Taylor's hypothesis and a relation between time 
and spatial spectra. Second order statistics are used to separate the deterministic kernel 
function and the weakly stationary driving noise. A non-parametric estimation method 
for the turbulence intermittency is suggested. 
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1. Introduction 

The Wold-Karhunen representation ([9], p. 588) states that every non-deterministic, 
one dimensional, stationary stochastic process Y = {Yh}heM.> whose two-sided power 
spectrum E(£) satisfyies the Paley- Wiener condition 

l0gEm dt<oo, (1) 



can be written as a causal moving average (CMA) 



Y h = (Y)+ / g(h - s)dW s , heR, (2) 

J — oo 

where W = {Wh}h&M is a process with uncorrelated and weakly stationary increments, 
(dWh) = and ((dWh) 2 ) = dh, with the angle brackets denoting the ensemble average. 
The kernel g is an element of the Hilbert space L 2 , i.e. \\g\\ 2 L 2 '■= J °° \g(s)\ 2 ds < oo and 
it is causal, i.e. it vanishes on (— oo,0). 

The converse is also true, i.e. every stationary stochastic process of form ^ satisfies 
0. If we consider h to be time, the representation ([2| has the physically amenable 
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feature of being causal, i.e. Y depends only on the past. The auto-covariance function 
of Y has a simple expression: for r G K, 



7y (r) = ((Y h+T (Y))(Y h (Y))) = / g(s + \r\) 9 (s)ds, (3) 

Jo 

while the two-sided power spectrum is 

E(0 = (2Mj Q oc lY ( s )co S (is)d S = \F{g(-)}\ 2 (0, where 



i r°° 

— / g(a)e-* a d8, t G K, 

'27T Jo 



denotes the Fourier operator. Condition (JlJ is crucial, since for a general stationary 
process a representation similar to ^ holds, but the kernel g may not vanish on (— oo, 0) 
and the integrals in ^ and ^ are extended to the whole real line ( 39 , Ch. 26). The 
rate of decay of the spectrum plays a pivotal role in determining whether a given process 
has representation ([2| or not, since excludes spectra decaying at infinity as exp(— £) 
or faster. In this context h or £ do not necessarily denote time or frequency. In the 
following x will denote a stream-wise spatial coordinate, K\ the associate wavenumber, t 
a time coordinate, and ui = 2nf the associated angular velocity with frequency /. 

In the turbulence literature the spectral properties of turbulent velocity fields were 
intensively investigated, starting from Kolmogorov's K41 theory [18-20 . In the Wold- 
Karhunen representation ^ the second-order properties of Y depend only on the function 
g, with no necessity to specify the driving noise W. This is analogous to K41 theory, 
where the second order properties of the velocity field can be handled without considering 
the intermittent behavior of the turbulent flow. 

From now on, we shall denote the mean flow velocity component by V. Moreover, 
we work with the usual Reynolds decomposition V = U + u, where U denotes the mean 
velocity and u is the time- varying part of V. Then Re := UL/v is the Reynolds number 
of the flow, by L we denote a typical length, and v is the kinematic viscosity of the flow. 

In K41 the first universality hypothesis claims that, for locally isotropic and fully 
developed turbulence (i.e. Re 3> 1), the spatial power spectrum El of the mean flow 
velocity fluctuations, has the universal form 

E l {kx) = {ev^f'HL^x) = v^LirjKi), (4) 

where r/ — (is 3 /s) 1 ^ 4 and v v — (ev) 1 ! 4 are, respectively, Kolmogorov's length and ve- 
locity, and < &l(-) is a universal, a-dimensional function of a-dimensional argument. As 
a cornerstone of the K41 theory, much effort has been devoted to verify Eq. Q and to 
determine the functional form of $l . For finite Re we can define $l in Q as the rescaled 
spectrum 

In an experimental setting with a probe in a fixed position, the spatial power spectrum 
can be estimated from the time spectrum E T L , using Taylor's hypothesis |35j . 

E L ( Kl ) = UE T L {U Kl ), (5) 

where ki = 2nf/U. This is regarded as a good approximation, when the turbulence 
intensity I = yftu 1 ) /U < 1 E2 EH]- Relation ((HJ) is a first order approximation, where 
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in general higher order corrections are feasible (T3] [53], and also error bounds can be 
obtained. 

Comparison of a large number of experimental low-intensity data sets |12l 1301 131] 
shows that in the dissipation range <&^ e is unvarying on a wide range of Reynolds num- 
bers. On the other hand, the inertial range does not exist for small Reynolds number [6], 
however, its length increases with the Reynolds number. The supposed infinite differen- 
tiability of the solution of the Navier-Stokes equation yields $^ to decrease faster than 
any power in the far dissipation range (nil] > 1) |37j . Moreover, the link between second 
order properties and third order ones is given by Kolmogorov's relation [TSl [2D 121] 

S 3 ( x ,t) = ~-e X + 6v— -- 4 J q z—dz, (6) 

where S n (x, t) = ((u(s + x, t) — u(s, t)) n ) is the n-th order longitudinal structure function. 
Formula (|6| has been used in [33] to determine an expression for the spatial spectrum 
which decreases like K"exp(— pK\) as K\ tends to infinity and a and j3 depend on the 
Reynolds number. Numerical [HI HH 1221 12H 132] and experimental [3D] |3T] studies con- 
firmed such exponential decay in the dissipation range. Under the hypothesis of local 
isotropy the local rate of energy dissipation is 

e x = lb V {d x V) 2 ^ l ^(d t V) 2 =e tl (7) 

where the approximation with the instantaneous rate of energy dissipation e t follows 
from Taylor's hypothesis. As customary in physics literature, we will drop the term rate 
for the sake of simplicity. The average energy dissipation can be calculated directly from 
the spatial spectrum [3] as 

e := (ex) = 151/ / 4E L ( Kl )d Kl ^ (e t ), (8) 

JR 

where the integrand k\El{k,i) for K\ > is the dissipation spectrum. 

The generality of the CMA model j2j requires in the context of turbulence modeling 
some interpretation of the parameters, in particular, the kernel function g. In this paper 
we estimate the parameters of model ([2]), bringing along some physical discussion to 
motivate our assumptions. 

In Section [2] the model ^ for the time- wise behavior of the time- varying component 
u of a turbulent velocity field is motivated first by an analysis of the literature yielding a 
discussion on the errors of Taylor's hypothesis and consequences for our model. Physical 
scaling properties of the kernel lead to a model for g in the inertial and energy containing 
range. Moreover, we suggest a deconvolution method to estimate the increments of the 
driving noise. In Section [3] the estimation methods of Section [2] are applied, using a non- 
parametric estimation of the kernel on 13 turbulent data sets, having Reynolds numbers 
Re spanning over 5 orders of magnitude. Moreover the increments of the driving noise 
for one of the considered datasets is recovered and analyzed. 

2. Time-wise turbulence model 

In this section we present the theoretical aspects of our turbulence model. Section 
2.1 is devoted to motivating the model by a discussion of Taylor's hypothesis and certain 
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refinements. In part Section 2.2 a rescaling of the kernel similar to Q is suggested, 
such that the rescaled kernel depends on the Reynolds number Re and the turbulence 
intensity 7 only. In part Section 2.3 we present a model for the inertial range and the 
energy containing range in order to deduce the scaling with Re of some features of the 
rescaled kernel. Finally, part Section 2.4 deals with the estimation of the increments of 
the driving noise. 

2.1. Does the Paley-Wiener condition hold? 

First note that an exponentially decaying power spectrum E violates ([!]) and, there- 
fore, it leads to a non-causal representation, regardless of any power-law prefactors. For 
spatial spectra this is not against intuition: the basic balance relations leading to the 
Navier-Stokes equation must hold in every spatial direction; on the other hand, the 
causality of representation ^ makes sense, when considering the time-wise turbulence 
behavior. 

It is known that the spatial spectra estimated via Taylor's hypothesis ([5| give larger 
errors in the dissipation range rather than in the inertial range, and that the error is in 
general positive, i.e. the time spectrum decays at a slower rate than the spatial spectrum. 
Lumley |23j derived from a specific model an ordinary differential equation relating time 
and space spectra. Lumley's ODE was solved in [6] and used on jet data with 7 = 0.30, 
resulting that the time spectrum obtained via ^ at KxT) ps 1 is 238% higher than the 
spectrum obtained with Lumley's model. Moreover, [13] showed that the power-law 
scaling in the inertial range is substantially left unchanged by Lumley's model. A similar 
effect has been already observed in |36j . comparing Eulerian and Lagrangian spectra; 
such spectra decay as w~ 2 and w~ 5 / 3 , respectively. 

Based on such facts we postulate that the time spectrum for every flow with turbu- 
lence intensity 7 > satisfies ([I]), and it is related to the spatial spectrum by 

E L ( Kl ) = UEi{ki(UK X )), (9) 

where A depends on 7, such that A/(ki)/ki — » 1 uniformly in K\ as I 1 0. The classical 
Taylor hypothesis in (|5| assumes that all eddies are convected at velocity U; however, it 
is likely that larger eddies propagate with a velocity of the order U, whilst the smaller 
eddies travel at lower velocity, resulting in a less steep decay [2 [35] . The function Aj 
accounts for such spectral distortion. In our framework we will not take the limit 710 
for two reasons: firstly, since the variance (it 2 ) is finite, 710 would mean that the mean 
velocity tends to infinity, and, secondly, if 7 1 0, in virtue of Taylor's hypothesis, the time 
spectrum (15| would decay exponentially like the spatial one, violating the Paley-Wiener 
condition pi). Therefore, in the considered setting, the limit 7 1 is singular, and we 
shall consider only the approximation for 7 <C 1. We shall do the same with the singular 
limit for Re tending to infinity ([TT], Section 5.2), denoting it by Re ^ 1. 

To date the resolution of experimental data is limited to scales in the order of rj but, 
if the data are not too noisy in the dissipation range, it is still possible to check, whether 
([!]) holds. In Figure [l^,) the time spectrum for the data set h3 is depicted; Figure [TJd) 
shows that the integral in ([!]) converges and that the dissipation range does not make 
any significant contribution to the integral. 
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2.2. Rescaling the model 

The CMA representation ^ for the time-wise behavior of the mean flow velocity 
component of a fully developed turbulent flow in the universal equilibrium range (i.e. 
inertial and dissipation ranges) can be rewritten as 



V t = U + C 2 [ gidit - s))dW S) te 

J — oo 



where g is a positive universal function, depending on the Reynolds number Re and the 
turbulence intensity /; C\ and C 2 are normalizing constants to be determined. The 
integral represents the time-varying part u in the Reynolds decomposition. The scaling 
property of the Fourier transform, i.e. that for c > 0, cT{g(c ■)}(£) = •7""{<7(")}(£/ c )s gi ves 

El(u J ) = (C 2 /C 1 ) 2 \F{g(-)}\*(u J /C 1 ). 

Using Q replacing $ L with i.e. without considering the limit for Re ^> 1, and (|9]), 
the relation between the rescaled spatial spectrum and the time spectrum is given by 

v^&t^m) = c/(c 2 M) 2 m9(-)}l 2 (A 7 (f/ Kl /Cx)) . (10) 




Matching both sides of (10 1, we get C\ = U/Qirrf) = f v , i.e. the Kolomogorov 
frequency, and Cf = U/v v e. The Reynolds decomposition can be rewritten as 

V t = U+J—[ g(f v (t-s))dW a , (11) 



where dW s = y/edW s such that ((dW s ) 2 ) = eds. The time-wise increments of the mean 
flow velocity, at time scale A, can be calculated from (11) as 

[g(t - s) - g(t - A - s)]dW s + / g(t-s)dW s , 

-00 Jt~A 

where g(-) := ^JU /v v g{f n )- If A -1 f t _ A g{t — s)dW s — > a.s. as A | 0, we have that 
the derivative process is, with the limit assumed to exist a.s., 

d t V = Urn A-HfV = J-f v f g'(f v (t - s))dW s , 
y v r) J— 00 

i.e. it is again a model of the form ([2| and we essentially exchanged integration and 
differentiation. Moreover, plugging Q into ^ we obtain 



15iA / s 2 $ L (s)ds, 



and, recalling the definition of v v and rj and (10), we get 



v 

ira-)}(A/(0) ■ \\J = 15-, (12) 



The Plancherel theorem, the fact that 



and (fl2|) give that Hs'll^ 2 ~ 15/2tt for I < 1. 
(I7| gives the instantaneous energy dissipation 

e ( := ^(9 t V) 2 = ^ jT jf* </(/„(* - s 2 ))g'{f n (t - Sl ))dW st dW S2 , (13) 

whose mean value is e = (et) = 15/27r||</||? 2 e ~ e when I <C 1. 

The constant f /ify is a-dimensional, serving as a scaling factor of the model. More- 



over, the rescaled model (111 indicates that {(dW s ) 2 } s£ r can account for the observed 
intermittency, i.e. it must provide all the higher order features of turbulence that can 
not be reproduced by a Gaussian model, as, for instance, the non-Gaussian behavior of 



the instantaneous energy dissipation {e t }t^M as indicated by (13 1. Moreover, we stress 
that g is a second order parameter of V, and so is g' . 

The model dW s = y/o^dB s with Brownian motion {i? s } se R has been suggested in 
[2], where yf<Js is the random intermittency process, assumed to be independent of B, 
and {cr s } s6 R is the instantaneous energy dissipation, with mean rate (a s ) = e. The 
major shortcoming of this model is that the Brownian motion assumption implies that 
the distribution of the increment process {S^V}teR is symmetric around zero for every 
scale A, which is against experimental and theoretical findings, especially at small scales 
(see e.g. [TT], Section 8.9.3). Such shortcoming amended in [3] by assuming the presence 
of a possibly non-stationary drift Z , in addiction to Y , where Z is smoother than Y (see 
e.g. [3J, Remark 6). The assumed smoothness of Z implies that the drift has a negligible 
effect on small scales increments, and is therefore more suitable for modelling phenomena 
at larger scales, such the energy containing range. 

From a second order point of view, everything depends on U, Re, e and u, which are 



the parameters of (11 1. Moreover, given that the dependence of g on Re is known, it is 
easy to simulate a process having the prescribed time spectrum. As long as only second 
order properties are of interest, one can indeed take {W s } s6 k = {£> s } s gk- If higher order 
properties are of interest, a realistic model for W is needed. 

2.3. Dependence of the kernel function on the Reynolds number 

It has been observed from experimental evidence (see e.g. [TT], Section 5.2) that the 
mean energy dissipation e is independent of the Reynolds number provided Re ^> 1. 

Since v^rj/v — 1, as in (4), U/v^ = Rerj/L holds. From Eq. (6.8) of [27| we have that 
U/v v cx Re 1/4 . From ((3), the variance of the mean flow velocity V is 

(u 2 ) = ((V-Uf)=v? ) \\g\\l 2 , (14) 

where v v is independent of Re, when Re 3> 1, and ||<7||l2 depends on Re (and to a lesser 



extent on /). Then the turbulence intensity of (111 is, using (14), 



r \fW) v, n _ 1/4 
1 = jj = -jjWWl* oc Re ' || 5 || L 2. 



G 



Since I must be independent of Re, we have \\g\\L 2 °^ Re 1 / 4 . 

A parametric model, suggested in [5] for the kernel g, is the gamma model 

g(t) = CRef*--^-***, t > 0, (15) 

where /iR C > 1/2 and (5r c , Cr > 0. This model yields a power-law time spectrum |27| 

B£( W ) - (27r)- 1 CLr 2 (MRe)(4 e + " 2 )- MRe 5 (16) 

where T(-) is Euler's gamma function. For instance, the von Karman spectrum [17] 
is a special case of (16), where /iR e = 5/6 and Cr e = \ArC/r(/iR e ) « 1.1431, and 



C is the Kolmogorov constant, which has been found to be around 0.53 over a large 
number of flows and Reynolds numbers, and therefore universal and independent of the 



Reynolds number |34j . For oj <C <5R e , the spectrum ( 16 ) is constant; hence u> rs <5R e can be 
interpreted as the transition frequency from the inertial range to the energy containing 
range; moreover, since the upper limit of the inertial range is independent of the Reynolds 
number, the size of the inertial range varies as J^ 1 . 

The gamma model has two essential shortcomings: firstly, it fails to model the steeper 
spectral decay in the dissipation range, but can be regarded as a good model for the in- 
ertial range. Secondly, the kernel has unbounded support, implying that the process is 
significantly autocorrelated even for large times, although the autocorrelation is expo- 
nentially decaying. 

From now on, we shall assume that the kernel function g has compact support, i.e. 
g(i) = for t > Tr 6 and t < 0, where Tr 6 is the decorrelation time [SB]. Since the inertial 
range increases with Re (see e.g. [27], p. 242) we expect i5r c to decrease. For the same 
reason we expect Tr to increase with Re. 

Explicit computations can be carried out for the truncated gamma model, considering 
a truncation at TR e and assuming that the failure of the gamma model in the dissipation 
range does not significantly affect ||<?||l2. Then 



\\g\\h =cL / Tr ° s^-Ve-^'ds = r (2mc = 1} = - 2 (2mc = 1; 2TrcSrc) (fee) 1 - 2 ^ 

\\y\\L Rey Q Cr 2 22airc -i v Ko ^ 



where T(a, z) := f°° s a 1 e s ds. If TR e <5R e 3> and it does not vary too much with 
Re, with the choice of the parameters as in the [yon Karman spectrum and the fact that 
\\g\\ L 2 cx Re 1 / 2 , we get S Rc cx Re" 3/4 . 

2.4- Filtering, driving noise and model building 

Let us now consider a second order stationary stochastic process Y with (Y) = 0. A 
filtered process X ([29 , Eq. 4.12.1) is defined as 



/•CO 

X t = / h{s)Y t ^ s ds 
Jo 
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where h is a causal function in L 1 . Invoking representation ^ and exchanging the 
integrals yields 



X t = h(s) / g(t -s- z)dW z ds 

Jo J -co 

ft f-t — Z 

I h{s)g(t — z — s)dsdW z 

1 -co JO 

(h*g)(t~z)dW z . (17) 



We shall assume that g belongs to the so-called Schwartz space S of rapidly-decaying 
function; i.e., smooth functions (j> : R M such that sup sgR \s p d q cf)(s)\ < oo for all 
integers p, q. The Fourier transform is well defined for all functions in S, it maps S 
onto S and the inverse Fourier transform of an element of S is again in S ([15]. Theorem 
7.1.5). Since in turbulence |^ r {<7( , )}| is assumed to decay faster than any exponential, it 
belongs to S, and consequently g € S as well. 



From ( 17) we know that the filtering affects the kernel function <?, but not the driving 



noise. Moreover, the spectrum of X is, using the convolution theorem, 

E X ( U ) = \H(g * ^)(-)}l 2 M = 2tt| 7{h(.)}7{ 9 (.)}| J H, 



If we know g, we can use (fl7|) to build a filter h such that X t sa S^W — W t — W t ^A- To 



this end we choose for appropriate a > 0, 

(h * g )(t) = ct> A (t) = c^ A e- |(t - A)trQ e(;)e(A - t), (18) 

Cq.a is a constant such that II^AOIlia = A, and O(t) is the Heaviside function. Then 
from Q follows that the filtered process X is a stationary process, with variance A and 
it is uncorrelated for lag r > A. Moreover, <^ A (i) tends to v A6o(t) as A J, 0, where So(t) 
is the Dirac delta in 0, justifying X « W t — W t -A as A <C 1. 



By the convolution theorem ( 18 ) is equivalent to 



F{Km<»)H9(-)}(o>) = HM-M"), (is) 

where saddle-point integration gives that |J r {0A(')}l( a - 1 ) asymptotically decreases like 
CJ -(a+2)/(2a+2) ex p(^— UJ a /( a + 1 ) ) as u qq anc j f or a fixed A > 0. We stress that the new 
kernel </>a(£) satisfies the Palcy- Wiener condition ([!]) for every a > and, therefore, it is 
a legitimate kernel function. 



W.l.o.g we assume that F{g (')}(<*>) ^ for all uj, we can deconvolve (19), to get the 
desired filter h. To ensure that h exist, it suffices to choose the parameter a such that 
F{<Pa(-)}(w) I J 7 {d( , )}( ul ) is in <5 and, as a consequence, that the inverse Fourier transform 
of such ratio exists. Moreover, if the ratio satisfies the Paley- Wiener condition ([T]), then 



h is causal. An ideal choice in the r.h.s. of (19) would be (h * g)(t) — 0(i)0(A — t), 
which can be retrieved for a I 0; however, since |.F{<7}|(u;) decays faster than any power 
of u, such h does not exist in S. 

The procedure shown in this Section is general, and it can be applied without restric- 
tion to any physically meaningful process, where a model of the form ^ applies. 
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hi 


2.72 


36.42 


0.76 


112 


5.74 


1.14 


20.56 


5.37 


0.76 


363.85 


0.98 


0.84 


h2 


0.52 


55.07 


1.11 


105 


2.87 


0.52 


19.19 


5.21 


0.73 


387.11 


1.05 


0.38 


h3 


20.14 


22.08 


0.55 


162 


22.97 


3.57 


21.47 


6.46 


0.85 


276.51 


1.38 


2.46 


h4 


0.85 


15.7 


0.49 


253 


5.74 


1.74 


24.04 


8.08 


0.81 


120.81 


3.06 


0.53 


h5 


13.14 


24.56 


0.66 


184 


45.94 


6.01 


10.98 


6.9 


0.95 


316.1 


1.25 


2.87 


h6 


10.98 


8.28 


0.36 


495 


22.97 


9.02 


23.34 


11.31 


0.83 


63.57 


5.95 


2.31 


h7 


64.37 


5.32 


0.26 


640 


91.88 


25.66 


22.58 


12.85 


0.87 


49.87 


7.53 


6.05 


h8 


1995.62 


2.25 


0.14 


930 


367.5 


175.68 


22.15 


15.5 


0.84 


22.95 


13.81 


37.37 


h9 


2120.81 


2.22 


0.14 


978 


367.5 


190.05 


21.64 


15.89 


0.84 


21.78 


15.79 


39.46 


hlO 


5192.01 


1.78 


0.11 


1005 


367.5 


285.13 


22.88 


16.11 


0.83 


18.13 


18.42 


60.46 


hll 


13580.72 


1.4 


0.1 


1336 


735 


570.02 


21.34 


18.58 


0.8 


11.76 


26.57 


108.72 


al 


1.38 


837.47 


140.01 


7216 


10 


1.21 


15.29 


43.16 


0.77 


0.55 


282.72 


1.97 


a2 


9.41 


442.42 


115.85 


17706 


5 


2.99 


28.23 


67.61 


0.81 


0.24 


882.39 


4.75 



Table 1: hl-hll: Helium data [7], al: Sils-Maria data [Ti| . a2: Brookhaven data \10\. The norm \\g\\ L 2 
is evaluated using (|14|. 



3. Estimation 

3.1. The data 

In this section we estimate the kernel function g non-parametrically, the parameters 



of model (11), and the increments of the driving noise W. Not surprisingly, the quality 
of the data in the present study does not allow us to perform a reliable estimation in 
the dissipation range. For the kernel in the universal equilibrium range, the gamma 
model ( |15[ ) is estimated by the non-parametric kernel estimation method as suggested 
in [5]. We analyzed 13 different data sets, whose characteristics are summarized in 
Table [T] Also dependence of the parameters on Re is considered. Two of the data sets 
examined here come from the atmospheric boundary layer (al P3] and a2 [TU]) and 
eleven from a gaseous helium jet flow (records hi to hll [7]). Since the data sets come 
from different experimental designs, characteristic quantities such as Taylor's microscale 
Reynolds number K\ = \J (u 2 )X/u based on Taylor's microscale A = \/l5(« 2 )^/e are 
considered. This Reynolds number is unambiguously defined and the rough estimate 
Ra oc \/2 Re holds ([27], P- 245). In all data sets the turbulence intensity / is rather high 
and, therefore, the expectation of the r.h.s of Q would not be a good approximation 
for e. Moreover, in most of the data sets the inertial range is hard to identify due to low 
Reynolds numbers. The mean energy dissipation e has been estimated from ([6]), following 
|21j . Estimation of e is a challenging and a central task, since the rescaling constants in 
( [TT| ) depend, in the K41 spirit, only on e, v and U. 

3.2. Kernel estimation 

The kernel function g characterizes the second order properties of the mean flow 
velocity V, i.e. spectrum, autocovariance and second order structure function. Classical 
time series methods have been employed to estimate the kernel function of high frequency 
data [5] with a non-parametric method. The method is essentially based on finding 
the unique function g solving ([3|, when the autocovariance function 7y is given. We 
estimate the discrete autocovariances 7y(«) := 7y(iA), where A is the sampling grid 
size and i are integer values. It has been shown [5] that the coefficients of a high order 
discrete-time moving average process fitted to such closely observed autocovariance, when 
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properly rescaled, estimate consistently the kernel function on the mid-point grid; i.e, 
g(i) « g((i + 1/2) A) for integers i. 

The results are presented in Figure [2^,) , where all the estimated kernels are plotted 
on a logarithmic time axis. Due to instrumental noise, the high turbulence intensity 
of the data sets and the non-infinitesimal dimension of the hot-wires, estimates in the 
dissipation range (tf n < 1) differ significantly from one data set to another. In general 
we notice that the steeper decay of the spectrum in the dissipation range is reflected 
by the fact that the kernel functions tend to for tf v -C 1, instead of exploding as the 
gamma model for /ir = 5/6. The only estimation of the kernel that is not perfectly 
aligned with the others is the one for h5, which looks slightly shifted to the left. This can 
be attributed to the difficulty in estimating e, and consequently f v . In Figure |2j;) the 
estimated \\g\\L 2 based on (14 1 are plotted against R\, and the power-law fitting shown 
in the figure returned \\g\\ L a = 0.5081 Ra 5 oc Re 25 . 

The decorrelation time Tfj e has been estimated by the first zero-crossing of the es- 
timated g. It increases with K\ (Figure [2^,)) and it follows empirically the law T Ko = 
0.1408 R A 13613 cx Re 6806 (Figure §;))• 

The transition frequency (5r g is obtained via least-squares fitting of the gamma model 



(15) to the non-parametric estimate of g for tf n > 5 (Figure |2p)). The statistical fits in 
Figure ^p) are very good for all data sets considered, at least when not too close to Tr c . 
The transition frequency follows the empirical law (5r c = 61.4917 Ra -1 ' 512 oc Re~ 0,7564 
(Figure [2p)), which is exceptionally close to the exponent —3/4 estimated for the gamma 
model, and <5R e ?Re is between 1.5 and 4.1 in all data sets. The estimated values of /^R, e 
are also close to the reference value of 5/6 « 0.834, with a mean value of 0.8218. The 
only notable outlier is the data set h5, which has been already proven to be somewhat 
anomalous. 

As said in Section 2.3, the truncated gamma model is not able to capture the sharp 
cutoff in the neighborhood of Tr c nor the rapid decrease in the dissipation range. To 
estimate how the variance of the model ([2]) is affect by the truncation at Tr c , we consider 
the quantity 

„ m v ||g(-)9(r Rc --)||| 2 r(2/3,2 Cl c 2 Re^) 

H(Re)= 1 rWaT 



where the latter is obtained using the gamma kernel ( |15[ ) with parameters //r c = 5/6, 
T Ro = cxRe 13 and ^Rc = c 2 Re Q , where C\,c% > and a > (3 > 0, as indicated by 
the least-squares in Figure [2j H represents the ratio between the variance ( 14 ) using 
a truncated gamma model and a non-truncated one, and it is a decreasing function of 
Re, and it tends to as Re 3> 1, indicating that the truncation is important when 
Re 3> 1; i.e., when T Re i5R e <C 1. Nonetheless, using the values of Ci,C2,oi,0 returned 
by the least-squares fitting, we obtain H w 0.991 for the dataset a2, which exhibits the 
highest Reynolds number among the considered datasets. Then we can conclude that, 
for the wide range of Reynolds numbers considered, the variance of the model ([2| with a 
kernel following the truncated gamma model does not differ in a sensible way from the 
variance of the classical von Karman model. 

Finally, the behavior of <5 Ro as function of Re estimated via the gamma model agrees 
significantly with the data, showing that the dissipation range does not contribute ap- 
preciably to ||<?||l2, which is in agreement with the idea that there is few energy in the 
dissipation range. 
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3.3. Noise extraction 



Relation ( 19 ) was calculated for a continuous time process with the goal of esti- 
mating discrete time increments of the driving noise. In practice, if we have data 
sampled with grid size A, we can reconstruct the sampled quantities in the frequency 
domain only on < / < F/2 = (2A) _1 . If A -C 1, we can use the approxima- 
tion J r {(/)A(-)}( W ) ~ 1 f° r w < tt/A; similarly, for the Fourier transforms the relation 
T{h(-)}(2TTk/(NA)) = ADFT{h A }(k) + o(A) holds fl2g], Eq. 13.9.6), for N/2 + 1 < 
k < N/2, h A = {/i(iA)}i=i jy, i.e., it is a function h sampled with grid size A and 
DFT is the discrete Fourier transform. The number of observations of the mean flow 
velocity V A is denoted by TV. Then the increments of the driving noise on a discrete grid 
can be recovered by computing 

where g is the non-parametrically estimated g, truncated at Tr c and properly zero- 
padded to increase its length to N. Moreover, IFT denotes the inverse discrete Fourier 
transform. For turbulence, the high-frequency condition A -C 1 can be regarded as 
satisfied whenever (2A) _1 = F/2 < f v j i.e., if the resolution of the data is of the order of 



the Kolmogorov's length. We apply (20) to the dataset h2, which shows low noise level 
and good resolution in the dissipation range (F/f v « 5.52); in Figure^) part of the 
obtained increments are plotted as an example. The increments of the driving noise show 
a clear intermittent behavior, exhibiting clustering and a non-Gaussian distribution. The 
clustering, i.e. the fact that larger increments are not isolated but appear in clusters, is 
reflected in the autocorrelation function of (S A W) 2 , showed in Figure |3ji), which is very 
similar to the one of e t (Figure |3j:)), defined in Q. It is also remarkable that those two 
autocorrelations are in the order of 10~ 3 for tf„ > Tr c = 105, suggesting that that V t 
may be independent of Vt- T for r > Tr c , rather than simply uncorrelated. Similar plots 
can be obtained for the other datasets as well. 

The similarity between those two autocorrelation function can be explained by the 
fact that the derivative of the kernel function, plotted in Figure [3^,), is concentrated in 
< i/„ < 2, it tends rapidly to zero for t J, and is small for tf„ > 30. Then we can 
heuristically think of approximate g'(t) with Q(t-f 1 )Q(f 2 - t) for some < Ti < T 2 -C 



Tr 6 , then (13) reduces to 



Moreover the increments of the intcrmittency process shows a density similar to the in- 
crements of the velocity, with positive skewness and exponential tails (Figure |4^) and d)) 
at small scales, and looking more and more Gaussian as the scales increase (Figure |4]>c) 
and Figure [4j3-f)). These phenomena are universally recognized as distinguished features 
of turbulence Ml. Ch. 8). 



4. Discussion and conclusion 

In this paper we proposed application of a large class of stochastic processes in the 
context of time-wise turbulence modelling. The class of CMA processes ^ is rather 
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flexible, with the only constraint of having a spectral density satisfying the Paley- Wiener 
condition ([!]) , which excludes processes with spectrum decaying too fast but it still allows 
them to have sample path infinitely differentiable with probability 1. Although causality 
is a reasonable feature of time-wise behavior of turbulence, it is not in the spacial one, 
and the link between the two is the often criticized Taylor's frozen field hypothesis. To 
the authors' knowledge, this is the first time the question is raised and it is worthy to be 
studied in more detail, e.g. using DNS simulation as in [T]. 

Essentially the CMA model distinguishes between second order properties, accounted 
by the kernel function, and higher order ones, depending on the noise, which can be spec- 
ified independently from each other, in agreement with the K41 theory. The dependence 
of the kernel function on the Reynolds number is analyzed, with a special regard to its 
behavior far away from the origin. The analysis in the time domain allows higher reso- 
lution of second order properties at higher lags, which corresponds to the leftmost part 
of the one-sided spectrum, showing in a clear way how the inertial range, proportional 
to <5rc, and the decorellation time Tr c increase with the Reynolds number. 

We propose a modification of the gamma model of [2] , with the parameters depending 
explicitly on the Reynolds number, modelling well the inertial range and the transition 
to the energy range. Unfortunately the present data, due instrumental noise, does not 
allow a precise analysis of the behavior of the kernel function near to the origin, which 
corresponds to the dissipation range. Such an analysis would be possible using data 
coming from computer simulations, and it is left for future research. 

Moreover, a method to recover the driving noise is proposed, essentially based on con- 
structing a linear operator to suppress the linear dependence of the data. The obtained 
noise, which is dimensionally the square root of the energy dissipation, shows some of 
the features of the energy dissipation e t , collectively known as intermittency. That shows 
that the second order dependence does not play an important role in determining the 
high order statistics. 

We conclude mentioning that the analysis performed in this paper holds in great gen- 
erality for one-dimensional processes. It is possible to model the full three-dimensional, 
time-wise behavior of turbulent velocity field with a similar CMA model, and under 
the hypothesis of isotropy a model for the three-dimensional kernel function can be ob- 
tained from a one-dimensional one, since the two point correlator depends only on the 
longitudinal autocorrelation function (see e.g. [27], p. 196). 
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Figure 1: a) Time spectrum for data set h3, the solid line indicates Kolmogorov's 5/3 law. b) Convergence 
of the Paley- Wiener integral | Jl} . The spectral density is estimated with the Welch method, using a 
Hamming window of 2 14 data points and 60% overlap. 
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Figure 2: a) Estimated kernels, against logarith mic time scale, b) Estimated kernels, plotted in log-log 
scale. Solid lines are the fitted gamma models ( | 1 5 p - c-e) Decorrelation times T^ e , transition frequency 
<5r 6 and ||g||^ 2 as a function of Taylor's microscale-based Reynolds number. Solid lines represents power- 
law fittings of the considered quantities vs. K\. 
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a) Derivative of the kernel function b) Estimated W t — W t -A 




Figure 3: a) Derivative of the kernel, b) Increments of the driving noise, estimated with formula | |20| . 
c) Autocorrelation function of the instantaneous energy dissipation e t . d) Autocorrelation function of 
the squared noise {Sf^W) 2 . The vertical solid lines in a), b), d) denote the decorellation time t = Tr c . 
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Figure 4: a-c). Histogram of the increments of the main-flow velocity for the dataset h2, respectively at 
scales I = A, 20A, 1000A. d-f) Histogram of the estimated increments of the intermittency process for 
the dataset h2, respectively at scales i = A, 20A, 1000A. 
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